Directed flow at mid-rapidity in event-by-event hydrodynamics 
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Fluctuations in the initial geometry of a nucleus-nucleus collision have been recently shown to 
result in a new type of directed flow (vi) that, unlike the usual directed flow, is also present at 
midrapidity. We compute this new vi versus transverse momentum and centrality for Au-Au col- 
lisions at RHIC using the hydrodynamic code NeXSPheRIO. We find that the event plane of vi is 
correlated with the angle of the initial dipole of the distribution, as predicted, though with a large 
dispersion. It is uncorrelated with the reaction plane. Our results are in excellent agreement with 
results inferred from STAR correlation data. 

PACS numbers: 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

Analyses of correlations between particles emitted in 
ultrarelativistic heavy-ion collisions at large relative ra- 
pidity reveal azimuthal structure that can be interpreted 
as solely due to collective flow Event-by-event hy- 
drodynamics 0-Q is a natural framework for studying 
collective flow. In each event, particles are emitted in- 
dependently according to some momentum distribution 
determined by a fluid freeze out surface. The most gen- 
eral azimuthal distribution can be written as a Fourier 
series with respect to the azimuthal angle (p of the parti- 
cle momentum: 
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where w„ is the magnitude of anisotropic flow 0, and 
is the reference angle in harmonic n. An equivalent 
definition is 



(2) 



where angular brackets denote an average with the prob- 
ability law dP/d(f). Once a convention is chosen for the 
sign of Vn, these equations define ^'n up to lixju. Both 
Vn and may depend on the transverse momentum 
and the rapidity y. 

Typically, the largest term in the series ([T]) is elliptic 
fiow, V2 @, and ^2 is usually referred to as the event 
plane. The only other harmonics measured so far at 
RHIC are v\ and W4 [oj. However, it has recently been 
pointed out that fiuctuations in the initial density pro- 
file generally result in different reference angles for 
every n. In particular, one expects additional triangular 
(^^3) and dipole (v\) [ll| components from fiuctua- 
tions, both uncorrelated with the event plane 4' 2. While 
triangular flow has already been studied in a hydrody- 
namic framework [5, 12, 13], there has not yet been a 
prediction for this new v\ that arises from fluctuations. 

The variation of directed flow with rapidity can be 



uniquely separated into even and odd parts: 
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where iif™" 



(y), Vl/evcn(y) and 

tions of w, while v\'^^{—y) = 
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are even func- 
Usual directed 



flow (9l.[T3-[T6| is vl'^'^ivi), and the corresponding is 
correlated with the reaction plane and 5*2. Teaney and 
Yan recently argued [ll| that fluctuations in the initial 
geometry create an even part which depends 

weakly on y. This even part does not contribute to ex- 
isting directed flow measurements, but can be isolated 
experimentally [T?| using an event-plane method where 
the weights are independent of rapidity (assuming that 
the detector is symmetric in rapidity). 4'i™"(?/) is pre- 
dicted to have little correlation with the reaction plane 
or 5*2, unlike It is also expected to have little 

dependence on pt or y. 

There is no dedicated analysis of wj'^''"(y) yet, but in- 
direct evidence has been obtained [l7[ from recent STAR 
correlation data [Tsj , and both the magnitude and pt de- 
pendence of v\ are in qualitative agreement with theo- 
retical expectations jTl| . 

We present the flrst quantitative predictions for tij™" 
in Au-Au collisions at the top RHIC energy, using the 
hydrodynamic code NeXSPheRIO ML Calculations of 



,,odd 



will be presented separately 



NeXSPheRIO 



solves the equations of relativistic ideal hydrodynamics 
using initial conditions provided by the event generator 
NeXus Fluctuations in initial conditions are stud- 

ied by generating many NeXus events, and solving the 
equations of ideal hydrodynamics independently for each 
event. NeXSPheRIO provides a good description of ra- 
pidity and transverse momentum spectra [2l|, and el- 
liptic flow vi [2^ . In addition, it reproduces the ridge 
observed in two-particle correlations [23], which is pro- 
duced by initial fluctuations, followed by collective flow. 
Teaney and Yan's dipole asymmetry is created by fluc- 
tuations followed by collective flow, much in the same 
way as the ridge [lO|, and should therefore be present in 
NeXSPheRIO. 
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NeXSPheRIO does not include effects of shear viscos- 
ity, which was recently implemented in event-by-event 
hydrodynamics [1| . While shear viscosity produces a siz- 
able reduction of elliptic flow [23, UHl , we expect its effect 
to be smaller on vi , following the general observation that 
damping is larger for higher harmonics [l3l |. 

In Sec. ini we explain how v°^'^'^ is calculated. In 
Sec. mil we present predictions for its pt dependence 
and centrality dependence in Au-Au collisions at the top 
RHIC energy and compare with existing STAR data. In 
Sec. IIV[ we study the correlation of directed flow with 
the initial dipole asymmetry defined by Teaney [ll[ . Our 
conclusions are presented in Sec. El We use vi, as a 
shorthand notation for vf"^'^, v[/cv<=n throughout this pa- 
per. 



II. MEASURING vi IN EVENT-BY-EVENT 
HYDRODYNAMICS 

The code NeXSPheRIO emits particles at the end of 
the hydrodynamical evolution using a Monte-Carlo gen- 
erator, and one can analyze events much in the same way 
as in an actual experiment. Particles are emitted inde- 
pendently. The only nonflow correlation in NeXSPhe- 
RIO is the correlation induced by resonance decays ^26*] . 
In principle, vi(pt,y) and ^i{pt,y) can be computed di- 
rectly from Eq. ^ by averaging over particles emitted 
from the Monte-Carlo generator. The only practical limi- 
tation is computer time, which results in finite statistical 
errors. In order to obtain reliable results, we integrate 
vi over the pseudorapidity interval — 1 < < 1, corre- 
sponding approximately to the acceptance of the time 
projection chamber of the STAR experiment 27[. We 



estimate by computing a weighted average ^ ep,i in- 
spired by the event-plane method used in experimental 
analyses [28| : 



Qcos^i EP.i = ^ Wj cos I 

Qsin'^ EPA = WjShu 

hl<i 



(4) 



where the sum is over particles emitted in the pseudora- 
pidity interval — 1 < 77 < 1, 0j are the azimuthal an gles 
of outgoing particles, Q > 0, and Wj is a weight [T^ 
depending on the particle transverse momentum pt : 



(5) 



where angular brackets denote an average over particles 
of all events in the centrality window and in the same 
pseudorapidity interval — 1 < 77 < 1. This particular 
choice of the weight eliminates both v°'^'^ in Eq. ([2]), 
and nonflow correlations from momentum conservation 
(though the latter are not implemented in this calcula- 
tion). We compute {pf) and (pt) by averaging over all 
events in the 0-40% centrality window, where pt spectra 



depend little on centrality [23]. w{pt) changes sign at 
Pt = {pl)l{pt) ^ 0.85 GeV/c. 

In order to achieve small statistical errors, we run the 
Monte-Carlo generator many times for each event. In 
practice, the summation in Eq. ([4]) runs over 7 x 10^ 
particles. This ensures almost perfect reconstruction of 
^'i for most events. We then measure 



Vi{pt) = (COS(0 - ^BKl)) 
Vla{pt) = (sin((?;> - *EP,i)) 



(6) 



where brackets denote an average over particles in a pt 
bin. If ^'i in Eq. ^ is independent of pt, then it coincides 
with '^EP,i, and vis{pt) vanishes identically. 
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FIG. 1. (Color online) Hadron v\{pt) and viaipt) in I77I < 1 for 
a single NeXSPheRIO event with impact parameter 6=1 fm. 
Results are shown with and without correcting for the net 
transverse momentum (see text). Vertical bars are statistical 
errors. 

Fig. [T] displays wi(pt) and wis(pt) for the first event 
we analyzed. vi(pt) has roughly the expected behavior: 
it vanishes linearly at low pt [301, changes sign around 
1 GeV [ill], and increases linearly at high pt [s?!. The 
solid line is a 3-parameter rational fit with the same prop- 
erties: 



vi{pt) 



aptipt - b) 



Pt 



(7) 



vis{pt) is not compatible with 0: it is in fact positive 
for all Pt , which means that there is a small nonzero net 
transverse momentum in the direction perpendicular to 
'ifEP,!- Averages in Eq. ([5]) are taken over the pseudora- 
pidity interval — 1 < < 1. 

This net transverse momentum could come from sev- 
eral sources. In principle, net transverse momentum can 
be generated by the initial dynamics, if a correspond- 
ing transverse momentum is transferred to "spectator" 
nucleons. An additional net transverse momentum can 
be generated in the pseudorapidity interval used here if, 
during the evolution of the system, net momentum is 
transferred to particles that end up outside this range of 
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pseudorapidity. In addition, however, there can be spu- 
rious momentum generated in NeXSPheRIO during the 
thermahzation stage, as an unphysical side affect of the 
transition from the microscopic model NeXuS to hydro- 
dynamics. A detailed explanation and the alternative 
procedure that exactly conserves conservation of energy 
and momentum is presented in Appendix \^ though its 
implementation and a thorough investigation of energy 
and momentum non-conservation in the standard scheme 
are left for future work. 

We are not aware of any measurement of the net trans- 
verse momentum of charged hadrons in a symmetric 
pseudorapidity interval. In particular, there is no exper- 
imental evidence that it is larger than statistical fluctu- 
ations. Since our vi results can be significantly modified 
by a net momentum, which is yet unknown experimen- 
tally, it is useful to present another set of results with a 
different (smaller) value of the net momentum. 

A simple way to correct for the net transverse momen- 
tum is by boosting the system in the transverse direction, 
and by adjusting the boost velocity so that the net mo- 
mentum vanishes in the new frame. This is difficult in 
practice, however, because the "system" here is the set 
of particles in the pseudorapidity interval — 1 < < 1, 
and 77 is not invariant under a transverse boost. We adopt 
the following approximate procedure: we define the boost 
velocity hy v = J2\ri\<iP/ J2\j]\<i ^- We then boost all 
the particles, select particles with < 1 in the new 
frame, and compute vi for these particles. Since this set 
of particles differs from the original set, the net trans- 
verse momentum is not strictly zero in the new frame, 
but it should be significantly smaller. Results after cor- 
rection are shown in Fig.[TJ uis(pt) is smaller in absolute 
magnitude, but still nonzero, because the correction is 
approximate. 

The next step is to average over events. This is done in 
the same way as would be done in an experiment. The 
scalar-product method [32| provides a simple prescrip- 
tion for doing the average, and it is known to give the 
same results as the event-plane method if the event-plane 
resolution is not too large Eq. (7) of [s^] gives: 



vi{EP} = 



(viQ) 



(8) 



where angular brackets denote average over particles in 
the numerator, and over events in the denominator. Due 
to the (almost perfect) event-plane resolution, there is no 
need to use subevents: the flow vector of one subevent is 
approximately half of the total flow vector, Qa ~ ~ 

g/2. 

Results presented in this paper are obtained from a set 
of 120 different NeXSPheRIO events, each corresponding 
to a different initial geometry. Each event is a Au-Au 
collision at 200 GeV per nucleon in the 0-60% centrality 
range. Events are uniformly distributed in centrality, and 
centrality is determined according to the number of par- 
ticipant nucleons. Compared with an actual experiment, 
we work with a limited number of events, but we generate 



a very large multiplicity in every event, so that our deter- 
mination of vi is not limited by statistical fluctuations. 
In practice, 75% of the events that we use have a res- 
olution parameter [sO] X = Ql larger than 3, 
implying an event plane resolution (cos A0) larger than 
0.97. Since x scales like Q, events with the lowest reso- 
lution have a negligible contribution to vi{EP}. 



III. PREDICTIONS FOR vi AT RHIC 
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FIG. 2. (Color online) wi(pt) and «is(pt) of hadrons in \ri\ < 1 
averaged over events in the — 40% centrality window. As in 
Fig. [1] results are shown with and without correcting for the 
net transverse momentum. Solid lines are fits using Eq. ([7]). 

Our results for vi{pt) and wis(pt) in the centrality bin 
— 40% are presented in Fig. [2l with and without cor- 
recting for the net transverse momentum in each event. 
After averaging over events, fis(pt) is much smaller in 
absolute magnitude than vi(pt). In fact, vis{pt) should 
be identically zero by parity symmetry. Any difference 
can be attributed to statistical uncertainty due to the 
finite number of events. Results in this centrality range 
are very similar with or without correction for the net 
transverse momentum. The pt dependence of vi is in 
qualitative agreement with theoretical expectations [ll| : 
it is linear and negative at low pt, reaches a minimum 
between —0.01 and —0.02, then increases and crosses 
at Pt ~ 1.1 GeV/c. 

Fig. [3] compares our hydro calculation for the central- 
ity range 20-60% to values of vi [17| inferred from STAR 
correlation data in the same centrality range .18j. As 
in Fig. [2I results are shown with and without correcting 
for the net transverse momentum in each event. For this 
centrality range, wi(pt) is slightly larger when the correc- 
tion is made, resulting in excellent agreement with data 
in the range < pt < 2 GeV/c. 

Fig. m displays the centrality dependence of — wi, aver- 
aged over hadrons with pt < I GeV/c, |ry| < 1. Results 
are smoother when NeXSPheRIO results are corrected 
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FIG. 3. (Color online) ui(pt) of hadrons in < 1 in Au-Au 
collisions (20-60% centrality) inferred from STAR correlation 
data [isl ] . Different sets of points correspond to different trig- 
ger particles and different assumptions concerning the «i of 
trigger particles (see O)- Solid curve: our hydro calculation 
averaged over events in the same centrality window, corrected 
for net transverse momentum (see text) and fit using Eq. ((TJ. 
The dashed curve represents the uncorrected value. 
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FIG. 4. (Color online) \vi \ — —vi of hadrons with < pt < 1, 
I?;! < 1 versus centrality. Dashed lines: raw NeXSPheRIO re- 
sults. Solid lines: after correcting for net transverse momen- 
tum (see text). Dotted lines: Rms dipole asymmetry (scaled 
down by a factor 10) versus centrality. Error bars on vi are 
obtained by assuming that the relative statistical errors on vi 
and £i™° are identical. 



CORRELATION WITH THE INITIAL 
DIPOLE ASYMMETRY 



Teaney and Yan [11 



have shown that fluctuations in 
the initial geometry of a nucleus-nucleus collision are ex- 
pected to create this new type of directed flow. Fluctu- 
ations break the symmetry of the initial density profile, 
and as a result there is, in general, one direction where 
the profile is steepest. This effect can be quantified as a 
dipole asymmetry in the initial density [ll| : 



eie 



i<E>i 



(9) 



where the averages in the right-hand side are taken over 
the initial transverse energy density profile, and (r, (p) 
is a polar coordinate system around the center of the 
distribution, defined by (re"^) — 0. If one chooses ei to 
be positive, then $i generally corresponds to the steepest 
direction for a smooth profile, and £i is the magnitude 
of the dipole asymmetry. In general, £i differs from — 
even at midrapidity — due to fiuctuations. 

For smooth initial conditions, one expects ^'i — $1 in 
each event. One also expects vi oc ei, in the same way 
as elliptic flow is proportional to the participant eccen- 
tricity [s^ l- Experimentally, vi will be extracted from a 
two-particle correlation, which scales like (tji)^, so that 
the experimentally measured vi should scale like the rms 
dipole asymmetry e™^ = ((ei)^) [S^- Fig. |4] displays 
the rms dipole asymmetry versus centrality. The val- 
ues are similar to those obtained using a Monte-Carlo 
Glauber calculation [11]. In particular, is in general 
larger for more peripheral collisions, and scales approx- 

— 1/2 

imately with the number of participants like A^p^i-t ; 
the same way as £3 [13] ■ The variation of — ui for pe- 
ripheral collisions is milder than that e™**. This can be 
attributed to the early freeze-out in peripheral collisions, 
which has an effect analogous to viscosity and breaks £\ 
scaling. 
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for the net transverse momentum. — ui increases mildly 
with centrality. The reason is that vi is created by fluc- 
tuations, which are larger for more peripheral collisions. 
The centrality dependence of vi is comparable to that of 




FIG. 5. (Color online) Distribution of the relative angle be- 
tween the event plane '^i and the initial dipole $1 for the 
centrality range — 60%. As in previous figures, results are 
shown before and after correcting for net transverse momen- 
tum. 
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With bumpy initial conditions, there is no simple di- 
rect correspondence between the initial geometry and the 
final momentum distribution Q ■ Fig. [5] displays the cor- 
relation between the angle of directed fiow, and the 
angle of the initial dipole asymmetry, $i. The correla- 
tion is quite clear, which confirms that the dipole asym- 
metry is a valid mechanism for generating vi. However, 
there is a large dispersion: some events develop directed 
flow in a direction very different from the initial dipole 
asymmetry. A similar dispersion has also been observed 
for elliptic flow between and $2 H and for trian- 
gular flow between 'i'^ and $3 36]. The dispersion of 
'I'l — $1 is qualitatively similar, but much larger. Note 
that in our calculation, Vfi is determined very accurately 
for each event, and this dispersion cannot be attributed 
to a finite event-plane resolution [28j . 

This large dispersion shows that, although the initial 
geometry specifies directed fiow and 'i'l completely, the 
information on directed fiow is not uniquely contained 
in initial dipole asymmetry. In the language of cumu- 
lants ^IL] , the dipole asymmetry represents only the first 
term in an expansion, and higher order terms contribute. 
We have also checked that the distribution of $1 is fiat: 
the direction of the dipole asymmetry is uncorrelated 
with the reaction plane (which is the x-axis in NeXSPhe- 
RIO), as expected [ll|. 
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FIG. 6. (Color online) Scatter plot of average vi in < pt < 
1 GeV/c, I77I < 1 and initial dipole asymmetry ei. Results for 
vi are corrected for net transverse momentum. 



Fig. m displays the values of fi, averaged over < pt < 
1 GeV/c, and the initial dipole asymmetry £1 for the 120 
hydro events used in our analysis. This figure shows that 
there is no one-to-one correspondence between vi and 
£1, though events with a large dipole asymmetry give on 
average a larger vi . 



V. CONCLUSIONS 

We have computed the directed fiow, vi, created by 
initial state fiuctuations near midrapidity using event-by- 
event ideal hydrodynamics for Au-Au collisions at RHIC. 
This Vi is an even function of (pseudo)rapidity and has a 
specific Pt dependence: it is negative below 1 GeV/c and 
positive above. Its minimum value lies between —0.02 
and —0.01 in the centrality range — 40 %. This is in 
contrast to the pt dependence of the usual, rapidity-odd 
directed fiow, which is typically iiegative for all pt at SPS 
(see the wi{3} resuhs in Ref. ^) and at RHIC 

The angle of directed fiow 5'i is correlated with the 
angle of the initial dipole asymmetry $1, but with a large 
dispersion. This confirms the idea of Teaney and Yan 
that the dipole asymmetry is the mechanism creating vi . 
However, it is only a rough picture, vi is not proportional 
to £1 on an event-by-event basis. More work is needed to 
understand how fiuctuations in the initial geometry are 
related to final momentum spectra. 

Our ideal hydrodynamics results are in remarkable 
agreement with preliminary experimental results [itI ITsj 
inferred from STAR dihadron correlation data in mid- 
central Au-Au collisions. This quantitative agreement 
suggests that vi is less sensitive to viscosity than V2 and 
W3. Detailed study of vi in viscous hydrodynamics is left 
to future work. Our results establish directed flow at 
midrapidity as a clear probe of hydrodynamic behavior. 
Experimental uncertainties on this vi are still large, but 
they can be significantly reduced by carrying out dedi- 
cated analyses at RHIC and LHC. 
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Appendix A: Initial conditions for event-by-event 
hydrodynamics 

In event-by-event hydrodynamics, initial conditions 
are taken from a microscopic model, which provides an 
energy- momentum tensor T^'^. Since local thermal equi- 
librium is usually not achieved in the microscopic model, 
rpp-u jg j^Q^ ^Yie energy-momentum tensor of an ideal fluid. 
By switching from the microscopic model to hydrody- 
namics, one must modify T**" — essentially approximat- 
ing a yet-unknown thermalization or isotropisation mech- 
anism. In doing so, one must choose quantities to keep 
continuous during this transformation, while the rest of 
the energy-momentum tensor necessarily changes discon- 
tinuously. 
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A common procedure is to transform T'^'^ into the lo- 
cal rest frame @ where the momentum density T*^* van- 
ishes. The fluid velocity is then defined as the velocity 
of the transformation, and the energy density is defined 
as e = T*^*^ in the local rest frame. Mathematically, this 
boils down to a diagonalization of [2] : is the nor- 
malized time-like eigenvector and e is the corresponding 
eigenvalue: 

Tf'^u" = eu^. (Al) 

The pressure P is then defined from e using the equation 
of state of the fluid, and the energy-momentum tensor of 
the fluid is defined as usual as = {e+P)ut'u'' -Pg>"' . 
Thus the local rest frame, defined by the fluid velocity 
u'^ , and the energy density in this frame are chosen to be 
continuous through the thermalization transformation. 

By transforming back to a common lab frame, however, 
it becomes apparent that the energy and momentum den- 
sity are changed in this process, in an uncontrolled way. 
For simplicity, consider the case where the transition be- 
tween the microscopic model and hydrodynamics is done 
at an initial time to in the lab frame. The total energy 
and momentum of the fluid are the integrals of T^" and 
T"* over the entire space. The above procedure changes 
T^" and T°*, and hence violates conservation of energy 
and momentum, even though the corresponding quanti- 
ties in the local rest frame T^^g^ = u^T^" are continuous. 

Alternatively, rather than demanding continuity of u^^ 
and e, one can instead choose a procedure that respects 
conservation of energy and momentum, which we de- 
scribe here. 

It should be noted that the above procedure is not al- 
ways used 3 , and in hydrodynamic calculations without 
a model for the initial microscopic dynamics (e.g., when 
only the transverse profile is postulated from a model, 
while the initial fluid velocity is set by hand) the issue 



does not arise d, HI . 

In general, the transition from the microscopic model 
to hydrodynamics is done across a space-like "freeze-in" 
hypersurface E. The total energy and momentum across 
E is 

P^= / T^"'da^, (A2) 

where dui, is the elementary time-like vector normal to 
E. 

Global conservation of energy and momentum de- 
mands continuity of P^ across the freeze-in surface. This 
can be ensured by demanding local energy/momentum 
conservation {d^T^'^ = 0), which requires continuity 
of the energy and momentum flux across the surface. 
Thus, to generate initial conditions for ideal hydrody- 
namics in a way that respects energy-momentum conser- 
vation laws one must demand local continuity of T^^'^da^. 
For example, in the case of a constant time surface, 
dui, oc (1, 0, 0, 0), this amounts to continuity of the energy 
and momentum density T'^^ Q. 

In general, this gives 4 continuity equations relating 
the hydrodynamic variables (e, u^, P) to the energy- 
momentum tensor T^'' from the microscopic model, 

T^'^dcr^ (e + P)u''u''da^ - Pg>'''da^. (A3) 

Further specifying an equation state P(e) uniquely deter- 
mines the hydrodynamic initial conditions, without any 
additional freedom. 

Note that this method agrees with the former method 
only at points along the frecze-in surface where da^ oc u'^, 
with it^ defined by Eq. (|Aip . Everywhere else along the 
surface, demanding conservation of energy and momen- 
tum results in values for e and in the hydrodynamic 
phase that are different from those acquired from the mi- 
croscopic model by Eq. (jAl[) . 
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